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Abstract We present a physical methodology to reconstruct the trajectory of 
interplanetary shocks using type II radio emission data. This technique calculates 
the shock trajectory assuming that the disturbance propagates as a blast wave in 
the interplanetary medium. We applied this Blast Wave Reconstruction (BWR) 
technique to analyze eight fast Earth-directed ICMEs/shocks associated with 
type II emissions. The technique deduces a shock trajectory that reproduces the 
type II frequency drifts, and calculates shock onset speed, shock transit time 
and shock speed at 1 AU. There were good agreements comparing the BWR 
results with the type II spectra, with data from coronagraph images, in situ 
measurements, and interplanetary scintillation (IPS) observations. Perturbations 
on the type II data affect the accuracy of the BWR technique. This methodology 
could be applied to track interplanetary shocks causing Til emissions in real¬ 
time, to predict the shock arrival time and shock speed at 1 AU. 

Keywords: Type II Radio Burst, Interplanetary Shock Waves, Coronal Mass 
Ejections 


1. Introduction 

Coronal mass ejections (CMEs) associated with interplanetary (IP) shocks are 
the most important phenomena of solar activity for space weather purposes. 
Together, IP counterparts of CMEs (ICMEs) and shocks are the main triggers 
of major geomagnetic storms when they interact with the Earth’s magnetosphere 
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{e.g. lOntiveros and Gonzalez-Esparza[ 120101) . Since fast Earth-directed ICMEs 
and associated shocks (from now on ICMEs/shocks) represent potential hazards 
to the geomagnetic stability, tracking them in the IP medium is a priority 
for space weather purposes. In general, the tracking of an ICME/shock event 
requires combining multiple observations, by different instruments, in order to 
approximate their trajectories in the inner heliosphere {e.g. IBisi et aO I20I0|I . 
However, despite the current observational capabilities, such tracking is still a 
very difficult task. 

The kinematics of ICMEs/shocks has a major interest in space weather stud¬ 
ies. It is well known that fast CMEs (> 900 drive shock waves near 

the Sun dVourlidas et ai, 2003D , and that IP shocks significantly decelerate as 
they propagate in the IP medium ( [Gosling et ai, 1968t . Furthermore, com¬ 
paring near-Sun (plane-of-sky) and in situ (at 1 AU) speeds, shows that fast 
CMEs decelerate as they propagate through IP space {e.g. [Lindsay et al.\ 119991 
[Gopalswamy et aH 120001 [Gopalswamy et 120051) . In general, the initial speed 
of fast CMEs as measured in coronagraph images is always faster than its 
corresponding transit speed at 1 AU, which in turn is faster than its in situ 
ICME speed at I AU. In fact, [Corona-Romero and Gonzalez-Esparza[ (I2011|) 
argue that this deceleration indicates that IP shocks are not driven any more by 
their corresponding ICMEs when they reach the Earth’s orbit. 

The speed of IP shocks also decreases with heliocentric distance ( Dryer[ 119741 


IPinterl[l9821 [Gopalswamy et a/.[ 120051 [Reiner, Kaiser, and Bougeret 120071) and 

shock’s speed apparently decreases with the inverse of the square-root of dis¬ 
tance ( [Dryer and Smart] 119841 ISmart and Sheal I1985|l . Such a decaying rule is 
a signature for blast waves in density profiles that decay with the square of 
distance, like IP medium does (see, [Rogers[ 119581 ICavaliere and Messinal I1970|l . 
Then, the deceleration of IP shocks was interpreted as a blast wave propa¬ 
gation {e.g., ISmart and Sheal II9851 [Pinter and Dryer[ II990|I which may start 
somewhere before I AU ( [Gorona- Romero and Gonzalez-Esparza, 201 1[ ). In fact, 
[Feng et ^(|20I0p empirically concluded that many of in sitw-detected (at 1 AU) 
IP shocks already evolve as blast waves. Furthermore, Pinter and Dryer] (|1990l) 


and [Corona-Romero, Gonzalez-Esparza, and Aguilar-Rodriguez (|20I3p indepen¬ 
dently found that trajectories of IP shocks are well described by blast wave 
propagation for heliocentric distances larger than ^ 25 — 45 i?©. Then, for the 
purposes of this work, we will assume that trajectories of IP shocks can be 
approximated by blast wave equation. 

Type II (TII) decametric-kilometric radio burst emissions is a useful signature 
to track the propagation of an IP shock ( [Pinter, 1982[ [Pinter and Dryer, I990[ ) . 
These events are characterized by a narrow band of intense radiation drifting 
to lower frequencies as time increases (distance from the Sun). TII emissions 
are produced by the excitation of plasma waves in the ambient solar wind 
by a disturbance (say an ICME/shock) propagating outward from the Sun 
( [Cane and Stone, I984[[Cane, Sheeley, and Howard, 1987 ). These TII emissions 
occur at the fundamental and/or harmonic of the plasma frequency (/), which 
is proportional to the square root of the local electron density (n), at the source 
region. 
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TII emissions usually start at frequencies below 150 MHz, where the distur¬ 
bance is just a few solar radii away from the Sun, and may extend down to the 
kilometric domain, slowly drifting to lower frequencies all the way to 1 AU, where 
the local plasma frequency of the solar wind is about 25 kHz. Type II emissions 
are classified according to their wavelength regime as: metric (m), decame- 
ter/hectometric (DH), and kilometric (km) bands. It is now well established that 
DH to km TII emissions are caused by the propagation of ICMEs/shocks through 
the IP medium | |Cane and Stone, I984[[Cane, Sheeley, and Howard, I987p . How-| 
ever, we should keep in mind that not necessarily all ICME/shock events generate 
Type H radio emissions dGopalswamy et ai, 2008D . 

Since beyond the supermagnetosonic point the solar wind density decays with 
the square of the heliocentric distance (n oc r“^), the frequency-drifting of TII 
emissions allows us to track the propagation of their source region. These radio 
observations can be used to reconstruct the propagation of an ICME/shock. 
IPinterl(II982p and |Pinter and Dryer] (| 1985)1 studied the propagation of IP shocks 
by using several techniques, including frequency drift of the associated TII 
emissions. Subsequently, [Pinter and Dryer] (|I990I) used data derived from TII 
emissions to develop a semi-empirical model to approximate IP shocks trajecto¬ 
ries and transit times from the Sun to Earth. In their model the shock trajectory 
has two stages: an initial short interval of constant speed, followed by a blast- 
wave propagation. Subsequently. IWatari and Detmanl(|1998p used this blast wave 
model to explore relations between transit times and arrival speeds (at 1 AU) 
of IP shocks. 

ILeblanc et al.\ (I200ip traced shocks/CMEs from the solar corona up to I 
AU, by combining white-light coronagraph, TII and in situ data. Afterwards, 
[Reiner, Kaiser, and Bougeret (|2007p proposed a kinematic model for the propa¬ 
gation of ICMEs/shocks which reproduces the TII frequency drift. This model as¬ 
sumes an arbitrary kinematics where initially the CME/shock suffers a strong de¬ 
celeration followed by a constant-speed propagation. [Cremades, St. Cyr, and Kaiser] 
()2007p suggested a method to predict the arrival of ICMEs/shocks associated 
with type H emissions assuming a constant velocity in the IP medium. On 
the other hand, following the work by [Reiner, Kaiser, and Bougeret] (|2007p . 


Conzalez-Esparza and Aguilar-Rodriguez] ()2009p suggested a technique to ap¬ 


proach the local shock speed using TII emissions. This technique assumes that 
the shock speed can be considered constant during short time intervals. By 
analyzing a few case events combining coronagraph, IPS, in situ and TII data, 
they found that fast shocks/ICMEs gradually decelerate from near the Sun to I 
AU. 


The aim of this paper is to describe a new physical reconstruction of ICME/ 
shock propagation using TII data. The methodology assumes that IP shocks 
decelerate as blast-waves and uses TII data to reconstruct the shock trajectory 
from near the Sun (~ QRq) to I AU. This reconstruction deduces: (I) the shock 
initial speed (which can be compared with the CME initial speed reported by 
SOHO/LASCO), (2) shock arrival speed and (3) shock transit time to 1 AU 
(which both can be compared with in situ measurements). We analyzed eight 
case events of fast-halo Earth-directed CMEs associated with kilometric TII 
emissions and ICME/shock counterparts. In general, we found a very good agree¬ 
ment comparing the results obtained by the blast wave reconstruction technique 
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and the observations. The outline of the paper is as follows: Section [5] explains 
the reconstruction technique using TII emissions; Section [3] presents the analysis 
of eight case events; and Sections [Hand [5] present the discussion and conclusions 
respectively. 


2. Physical reconstruction of shock trajectories using Type II data 

The methodology to obtain the trajectory of ICME/shock events associated 
with TII emissions can be summarized in three steps as follows: (1) from the 
TII data we obtain the central emission frequencies and their associated band- 
widths; (2) we apply a blast wave reconstruction model to the shock propagation 
using TII data; (3) we compare directly the shock parameters obtained by the 
reconstruction with chronagraph and in situ measurements. 

2.1. Type II central emission frequency 

In practice, most observations of TII radio emissions are complex and show 
important fluctuations, and even gaps, in intensity, frequency emission and drift 
rate. Most of these fluctuations derive from the solar wind structure, the IP 
magnetic field topology and the shock front evolution ( [Knock and Cairns, 2005[ ) . 

Due to such complexity, one of the problems in analysing TII emission is to 
identify the central emission frequency ((/)) of the TII radio burst and the asso¬ 
ciated bandwidth (A/). In order to simplify this task, we apply the technique by 
[Aguilar-Rodriguez et a/.|(|2005|l and jGonzalez-Esparza and Aguilar-Rodriguez|(|2009|l | 
to infer (/) and its associated A/. We present in Figured] a case study to illus¬ 
trate this technique. Figured^ shows the radio spectrum of event 1 (see Tabled]) 
detected by Radio Receiver Band 1 (RADI) http://lep694.gsfc.nasa.gov/waves/data_products.htmLI 
RADI is one of three radio receivers of the WAVES experiment on board the 
WIND spacecraft ( ]Bougeret et al., 1995] ). RADI covers a frequency range from 
20 to 1,040 kHz in a frequency band divided into 32 real channels. The frequency 
emission is observed at 15:00 UT on 6 June at 600 kHz, drifting slowly down 
to 60 kHz at 00:30 UT on 8 June. The TII emission is characterized by a single 
tone with intervals of very intense emission, which lasts several hours. 

In Figure db we isolated the TII event from the dynamic spectrum by setting 
to zero any emission outside the TII emission feature. Figure dt* also shows a 
vertical dashed-line indicating a time sample (04:24 UT) on the TII emission. 

Figure dt shows this sample as a flux density vs frequency plot (asterisks), it 
also presents the Gaussian fit (dashed-line) to the flux-density distribution. This 
Gaussian approximation computes the values for (f)i (dash-dotted line) and 
Afil2 (dotted-lines) at the given time (U) (where i is the index in the data 
series). The systematic repetition of this process, at every time sample of the 
spectrum, gives us also criteria to determine whether the time sample in the 
spectrum has a well-defined central emission or not {e.g., there are some spectra 
with multiple maxima in the flux density distribution, for these cases we discard 
the data point). After applying this filtering technique to the TII data, we obtain 
a set of m values for {f)i, A/i/2 and U, with 1 < i < m. 
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Figure 1. Case study, central frequency and associated bandwidth deduction for event 1 
(LASCO/CME on 6 June 2000). (a) WAVES/RADI dynamic spectrum associated with the 
propagation of the ICME/shock. (b) Isolated TII emission between 6-Jun-2000 22:00 UT 
to 7-Jun-2000 06:00 UT. (c) Flux density vs frequency for 1-min RADI spectrum at 04:24 
UT. White dashed-line in panel (b) marks a temporal sample (at 04:24 UT) to perform the 
frequency analysis and the dashed black-line, in panel (c) points out the gaussian fit. Further 
description of the figure in the text. 


ICorona-Romero, Gonzalez-Esparza, and Aguilar-Rodriguez] ((2()13|) simulated 
TII emissions using a physical analytical model of ICMEs/shocks propagation. 
These synthetic radio bursts did not follow the central frequency of the associated 
TII emission but their lower frequency edge (/). Assuming that the shock front 
is the TII source emission, the shock leading front is at the farthest distance from 
the Sun, which is associated with lower solar wind densities emitting at lower 
frequencies. Based on this previous experience, we will use the lower frequency 
edge (/) of the TII emission for the regression model. From the analysis of radio 
data commented on before (see Figured]), we define the lower frequency edge as 
/■ = {f)i ~ As we will explain below, with these frequency data we can 

reconstruct the trajectory of the IP shock. 
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2.2. Blast Wave Regression technique 

The Blast Wave Regression (BWR) technique reconstructs the trajectory of 
shock fronts using Til data. In order to do so, this technique combines Til 
frequency data with the kinematic equations of a blast wave. A blast wave is a 
shock wave which has a hnite amount of energy, the shock is transferring energy 
to the medium, decaying continuously as it propagates. In the case of an IP 
shock, a blast wave propagation implies that the shock decelerates following a 
particular rate ( [Rogers, 1957p . On the other hand, a driven shock is a shock 
wave that has a continuous source of energy from a driver, this source input 
energy replaces the energy that the shock transfers to the medium. In the 
case of an IP shock, a driven shock propagation implies that the shock and 
its driver (ICME) propagate at almost at almost the same speed ( [Parker, 1961D . 
Previous studies, using in situ data, suggest that IP shocks are not driven by 
their ICMEs further away from 1 AU ( [Burlaga et al., 1981[[Feng et al., 2010] ) . In 
fact, [Corona-Romero and Gonzalez-Esparza (|2011l) claim that the CME/shock 
driven stage lasts just a few hours from which the IP shock propagation follows 
a blast wave decay. We should keep in mind, as commented on before, that the 
deceleration of fast CMEs/shocks in the IP medium is well established. For this 
study, we assume that IP shock waves follow a blast wave propagation from 
a certain heliocentric distance, which on average is about 15 Rq (this value is 
obtained from the reconstruction algorithm). 

ICavaliere and Messinal (|1976|l investigated the self-similar solutions for shock 
waves propagating through different density profiles. According to them, when 
density decays as r~^, the speed of a fast blast wave decreases as 
with T a constant of proportionality that modulates the decreasing rate of the 
shock speed. Furthermore, [Pinter and Dryer[ ([199011 found that shock speed is 
convected by solar wind expansion. Thus, the speed (v) of an IP blast wave can 
be approximated by: 



V = {vo- Wsw 


( 1 ) 


T 


where vq is the initial shock speed (at t = t) and Wsw the solar wind speed. We 
can express the leading front position of a blast wave (r) by integrating Equation 
© over time: 



with To being an integration constant related with the initial position of the 
shock front (at t = 0). Notice that Equations ([T]) and (HI are valid for t > t. 
In order to simplify the integration of Equation we assumed v = vq for 
0 < t < r. 

Then, from the relation between plasma frequency, solar wind density and 
heliocentric distance, we can define the position of source region by: 


n = r{U) = QNy/niAUY ■ 


( 3 ) 
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Where uiau is the solar wind density at 1 AU and Q = 

with e, £o, and rUe being the fundamental charge, electric constant and the 
electron mass, respectively. Additionally, N indicates the fundamental frequency 
(= 1 ) or the first harmonic (= 2 ) and / is the lower frequency edge. 

Now, by assuming that the location of TII source region is near the shock 
leading edge, we combine Equations © and to get: 




2/3 

^ 3^2 H " ^2^1 ^0 


(4) 


Equation o is a dispersion relation between frequencies and times that allows 
to reconstruct the kinematic evolution of IP shocks. The constants 03 , 02 and 
ao are related with the blast wave equations by: 


'^sw — 

r = 

^0 = 

We calculate the values of 03 , 02 and oq employing a gradient-expansion regres¬ 
sion algorithm ( |Bevington and Robinson, 2003[ ). Using the regression algorithm 
we can describe the kinematics of shocks by Equations © and as well as 
simulate the radio TII emission by combining Equations ([3]) and ([2]). In the 
Appendix we present the relations to estimate uncertainties associated to the 
BWR technique. 


asy/niAuQ N, 


asro - aoWsw 


3/2 


Ws 


^2^sw 


1 _L - 

1 + o - 

803 


(5) 

( 6 ) 

(7) 


3. Case events 


We selected eight case events from the WIND/WAVES type II bursts and CMEs 
catalogue (http://cdaw.gsfc.nasa.gov/CMEJist/radio/waves_type2.html) during the| 
decaying phase of solar Cycle 23 (2000-2005). Table[T]shows the event list. In all 
cases, the associated TII emissions required to be clearly identified and to last 
several hours. All events were catalogued as halo CMEs (Earth-directed), with 
(plane-of-sky) initial speeds above 900 kms~^ (Tabled! column 3), and were asso¬ 
ciated with solar flares, as reported by SOHO LASCO CME dGopalswamy et al, 2009D .| 


The flares’ locations were near the Sun’s center (see last column of Table dj) to 
reduce possible directivity issues between the observer and TII source region. All 
the shocks and ICMEs were reported in IR ichardson and Canel (|2()1()|1 catalogue, 
but event 4 which possibly was missed due to a data gap in ACE registers 
(see IJian et all 120061 catalog). We identified event 4 from WIND in situ data. 
Column 6 in Table [1] shows the transit times (from near the Sun to I AU) of 
shocks and ICMEs. Column 7 shows the speeds of shocks and ICMEs at I AU 
obtained from WIND data. We estimated the shock local speed at I AU {viau), 
applying velocity coplanarity on the upstream and downstream in situ data. We 
calculated the ICME speed averaging WIND in situ data on the ejecta leading 
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edge. Since there was a data gap in event 5 impeding the use of the velocity 
coplanarity to estimate the shock local speed, for this case, we used instead as 
a proxy the shock transit speed to 1 AU. 

The first step in the methodology was to study 1-min resolution TII frequency 
data following the procedure described in Figure[I] Subsequently, we reduced the 
frequency data temporal resolution to 0.5 hours, through an averaging process, 
in order to smooth out frequency irregularities to facilitate the convergence 
of the reconstruction algorithm. The averaging frequency data is presented in 
the bottom panels of Figures [2] and |3] as open-orange squares. Afterwards, we 
proceeded with the BWR technique by setting ro = 6 Rq (LASCO C2 external 
edge) for all the case events. The value of uiau (see TablelU column 4) was taken 
from in situ measurements 8-12 hours before the shock arrival. The initially 
supplied function for the regression algorithm used vq = Vcmeo, Wgw = Wswi 
and T = 5 ro /VcmeO (see Table [T|) . Additionally, we set A^ = 2 for all the events 
(hrst harmonic frequency) and we used the detection time of each event as t = 0 
(TablelU column 2). 

3.1. Results 

We present the results of the case events in Figures [5] and [3] In both figures, 
each column represents a single event, meanwhile the rows, from top to bot¬ 
tom, show the speed and heliocentric distance of the shock fronts in time and 
the dynamic radio spectra. The vertical dash-dotted lines point out the event 
detection-time reported by LASCO CME catalogue (see Tabled! column 1). 
Upper and middle panels of the figures present the shock trajectories (speed 
and position versus time) obtained from the BWR technique (solid-red lines). 
In the panels, we point out the initial (fill circles) and final (open circles) 
parameters calculated by the BWR technique. In these panels, for compar¬ 
ison with the BWR results, we also present data from: coronagraph images 
(open-diamonds) as reported by SOHO LASCO CME ( jCopalswamy et ai, 2009[ ) 
catalog; in situ data of ICME (down-ward open-triangles) and shock (asterisks) 
arrivals at 1 AU; and IPS speed measurements (open-stars) as reported by 
Gonzalez-Esparza and Aguilar-Rodriguez| (120091) and IBisi et aU (|2010ll . In the 
bottom panels of the hgures, we overplot the simulated radio emission (solid- 
white line) obtained from the BWR technique as well the radio data used in our 
analysis (open orange-squares). 

In most cases, the simulated radio emission approximates very well the fre¬ 
quency data. For events 6 and 7 (Figure|3|) the limited TII data derived into poor 
BWR results. However, it is remarkable that, despite the limitation in the TII 
data, it was possible to apply the BWR technique in both events. Furthermore, 
in all cases the simulated frequency drift systemically limit the lower edge of the 
observed TII radio emission. 

In general, the resulting trajectories (position and speeds) are in good agree¬ 
ment with coronagraph, IPS and in situ measurements. Table [3] summarizes the 
main results of the analysis of the 8 case events, comparing the BWR calcula¬ 
tions with coronograph and in situ data (Table I). According to Equations ([3]) 
and (HD, the calculated trajectory of the shock starts at t = r, time at which 
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Flare 

Location^-Class-^ 

1 

20000606-15:54 

1119 

4.7 

540 

41.2/44.1 

871/770 

N20E18-X2.3 

2 

20000714-10:54 

1674 

3.5 

700 

27.5/31.9 

1120/1040 

N22W07-X5.7 

3 

20010426-12:30 

1006 

2.5 

440 

40.5/49.5 

820/730 

N20W05-M1.5 

4 

20010924-10:31 
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3.2 
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34.5/39.0 

641/530'’ 

S16E23-X2.6 

5 

20011104-16:35 
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7.5 
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33.5/34.9 

31206/750 

N06W18-X1.0 

6 

20011122-23:30 

1437 

5.5 

419 

31.5/38.5 

1108/1040 

S17W36-M9.9 

7 

20040725-15:06 

1333 

1.8 
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31.9/35.1 

1153/1000 

N08W33-M1.1 

8 

20050513-17:12 

1689 

3.5 

410 

33.5/36.9 

'*1111/950 

N12E11-M8.0 


“ SOHO LASCO CME catalog ( http;//cdaw.gsfc.nasa.gov/CME_list/}. 
^ Detected in situ by Wind (http://omniweb.gsfc.nasa.gov/}. 

^ Calculated by velocity coplanarity. 


^ Taken from [Richardson and Canel II201QII catalogue. 

® SolarMonitor.org (http://www.solarmonitor.org/}. 

^ NO A A classification. 

^ Estimated value viau = (1 AU — ro)/TTsh due data gap. 
^ Value taken from f^si et a/.ll|2QlQ[l due data gap. 


CD 


Table 1. List of events (chronograph and in sztu measurements). Prom left to right: event number, date-time of CME 
detection, CME initial speed from LASCO observations; in situ values for proton density and bulk speed of solar wind; 
shock (/ICME) transit times and shock (/ICME) local speed at 1 AU; and location (-class) of the associated solar flare. 
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Table 2. Results from BWR analysis. From left to right: event number; values of r and associated! 
heliocentric distance calculated value of the ambient solar wind {wsw), calculated value of the 

shock initial speed (uq); calculated shock transit time (TT^^); calculated shock local speed at 1 AU 
{viAu)'i proportional absolute difference between calculated and in situ registered values of TTsh\ 
and viAu'i angle (a) between Sun-Earth line and location of associated solar flare. Column 
also shows the direct difference between calculated and in situ registered values of TTgh- 
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begins the shock blast wave propagation (Column 2). In general, the calculated 
shock initial positions {t ^ r) were located within the coronagraph (LASC0-C3) 
field of view. The calculated shock transit times (column 5) were quantitatively 
consistent with the in situ arrivals of the IP shocks. Column 7 lists the relative 
differences between calculated and in situ measurement of shock transit times 
(TT), where the average difference is about 10%. Column 8 lists the relative 
differences between calculated and in situ measurement of shock speed at 1 AU 
{viAu)- The average relative normalized difference was about 14%. Note that the 
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difference for event 4 is the largest; however, in this case, the estimation of the 
in situ shock speed has a large uncertainty. The fourth row of Figure [2] (upper 
panel) shows the largest error in the shock speed estimation. Neglecting event 
4, the average associated difference of shock local speeds drops to 9% (similar 
to the one for TTsh)- 

The BWR calculations of the shock speed evolution in Figures [2] and [3] (upper 
panels) present an initial strong deceleration that rapidly decreases with time. 
Such deceleration considerably reduces (^ 50%) the speed of IP shocks when 
they arrive to 1 AU (see Tabled columns 3 and 5). However, note that the 
deceleration is different for each event. We can estimate the decaying rate of 
shock speed by Equation m- 


dv oo / t \ 

where oq = [uq — Wsw\/t can be interpreted as an effective acceleration. Then, 
larger values of uq, imply stronger decelerations acting on blast waves. Events 2 
and 8 exemplify the effects of different values of vq — Wsw with similar r. For event 
2 (oo = 699 ms~‘^) the decrease in speed at 1 AU was ~ 50%; whereas for event 8 
(tto = 932 ms“^) it was ~ 61%. On the other hand, events 5 and 6 show the case 
for different values of r and close values of vq — Wsw For event 5 (oq = 81 ms~^) 
the shock speed decreases ^ 33%; whereas in event 6 (oq = 113 ms“^) the 
reduction was ~ 42%. 

The BWR calculated shock parameters at 1 AU for events 3,5,6 and 7 present 
differences with in situ data. Such differences may derive from contamination of 
the Til radio spectrum and/or directivity issues. On one hand, events 3 and 5 oc¬ 
curred during high solar activity intervals with occurrence of other CMEs before 
and after. This could lead into multiple (complex) CMEs propagating through an 
unstable ambient solar wind. This scenario might contaminate the TII dynamic 
radio spectra, misleading the results. On the other hand, the associated-flare 
locations of events 6 and 7 have a large angular distance (^ 35°) with respect to 
the Sun-Earth line (see Table [2j column 8). Using the flare location as a proxy 
for the CME propagation direction, we can estimate if the source emission is 
propagating far away from the Sun-Earth line. For these cases, the radio source 
may not be directly seen by the observer. TII emissions are affected by the line of 
sight between the observer and the source ( [Knock and Cairns, 2005P . This might 
be the case in event 7, where the TII emission does not appear as a continuous 
frequency drift, but like a chain of consecutive bursts. For events 6 and 7, this 
directivity issue might explain the reduced number of TII radio data, resulting 
into poor BWR results. In Figure [3] the simulated TII emission of events 6 and 
7 barely follows the TII spectral data. 


4. Discussion 

We present the BWR technique to reconstruct the evolution of shock fronts 
using TII emission data. The technique is based on three main assumptions: (I) 
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solar wind density decays as r~^ and expands with constant speed (wsw)', ( 2 ) 
the Til emission source is located somewhere near the leading edge of the shock 
front; and (3) IP shocks evolve as blast waves. According to the first assumption, 
this technique would be more suitable for quiet periods of solar activity, when 
isolated fast CME events propagate through a quasi-stationary ambient solar 
wind. On the other hand, note that the average initial positions of blast wave 
propagation (r,-), given by the BWR results, were located around 21 Rq (see 
Table [5] column 2) which is beyond the supermagnetosonic point. 

With respect to the second assumption, as some authors suggest, the source 
region of Til emissions might be located far from the shock’s leading edge 
( [Knock et al., 2003| ). This source location may affect the results of this tech¬ 
nique. The two events in the list whose associated flare locations were more 
than 30° from the Sun center (events 6 and 7 in the last column of Table 2) had 
significant differences with respect to the in situ measurements. It is possible 
that these errors were related to the location of the source regions, and that the 
BWR technique is more suitable for Earth-directed events, where source regions 
are close to the Sun-Earth line. 

It is not clear where ICMEs get exhausted of driving IP shocks in the inner 
heliosphere. However, beyond that critical distance (r,-), the IP shocks should 
evolve as blast waves, as the third assumption dictates. As commented in the 
Introduction, there are multiple works that suggest that IP shocks evolve as blast 
waves. For example, [Burlaga et al.j (119811) . by analyzing ICMEs/shocks detected 
by Voyager 1 and 2, found that the shock might be no longer driven by ICMEs 
at heliocentric distances around 2 AU. Departing from geometrical analysis, 
[Feng et aL\ l|2010|) found that at least 34% of shocks detected in situ (at 1 AU) 
were not driven by their associated ICMEs. [Corona-Romero and Gonzalez-Esparza] 
(|2011|1 studied the momentum fluxes between ICMEs and plasma sheaths using 
numerical simulations and concluded that IP shocks are no longer driven by 
ICMEs when they reach 1 AU. The results from our study suggest that IP shocks 
evolve as blast waves beyond ~ 21 Rq (t > t). This value is consistent with the 
ones separately found by Pinter and Dryerj (I1990|l (~ 0.12 AU = 25.8 i?©) and 
[Corona-Romero, Gonzalez-Esparza, and Aguilar-Rodriguezj (r2013p (^ 44.7 i?©) 
through empirical and analytic approaches, respectively. Furthermore, by as¬ 
suming blast wave propagation for IP shocks, the BWR technique was able 
to reconstruct speed and position profiles that showed qualitative and quan¬ 
titative agreements with a number of different data. Furthermore, blast wave 
propagation explains the well known fact of deceleration of IP shocks, as well as 
the unexpected growth in stand-off distances found by [Maloney a nd Gallagher[ 

iinni). 

In Figures 2 and 3 (upper and middle panels) we plot the shock trajectories 
using the approximation by [Reiner, Kaiser, and Bougeret[ ([200711 (dotted grey 
line). They proposed a two stages kinematic model to proxy the trajectory of fast 
ICMEs/shocks. Such a model consists in a pristine period of intense deceleration 
followed by a constant speed stage. Then they search for the trajectories of 
ICME/shocks that simultaneously satished data from coronagraph images, in 
situ registers and TII emissions. Comparing our solutions with the ones obtained 
by [Reiner, Kaiser, and Bougeret[ ()2()()7I1 . we notice similarities between them. 
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Blast wave speeds present an initial strong deceleration which rapidly decreases 
with time. This behavior resembles the one presented by the two stages kinematic 
model. Event 1 is an example where there is a good agreement between the two 
methodologies which both predict very similar shock arrival times to 1 AU. 

Though BWR technique is presented as a tracking method for IP shocks; it 
might be useful for forecasting proposes. According to our results, BWR is more 
suitable, but not restricted, for Earth-directed shocks occurring during quiet 
periods of solar activity. Our calculated shock trajectories showed a tendency 
to arrive later than their in situ counterparts. The average difference between 
calculated and in sztw-registered transit times was 3.5 h, a value similar to 
others 7h) derived from empirical (e.ff. |Cremades, St. Cyr, and Kaiser] 120071 
Feng et ^120091 IMostl et fl71l201i1l and analytical {e.g. Feng and Zhao [20061 


Songl [2(71(1 arrival-predicting tools. An important limitation of BWR technique. 


for predicting proposes, is the acquisition of in situ data. This data is taken 
around 8-12 hours before shock arrival, leaving few hours to issue alarms. Though 
this limitation can be easily solved by using IPS measurements or average values 
of solar wind density at 1 AU, exploring this topic escapes the scopes of this work. 

The BWR calculated shock speed profiles were consistent with coronagraph, 
in situ and IPS data. However, although our results showed good agreement 
with different data, we have to point out that those data correspond to different 
but related structures. On one hand, coronagraph images give information about 
(plane-of-sky) propagation of CMEs in the solar corona. On the other hand, IPS 
speeds might be related to plasma sheaths behind the IP shocks, and the in situ 
measurements to local characteristics of ICMEs and shocks at 1 AU. 

The BWR technique is sensitive to contamination and directivity issues on Til 
emission data. The applicability of BWR, and the certainty on its results, will 
depend on the availability and quality of the TII data. However, the calculated 
transit times and arrival speeds of shocks were quantitatively consistent with 
their in situ registered counterparts. In average, our results showed ~ 90% of 
agreement with in situ data (see Table [5J column 7). Even for unfavorable events 
(say events 3, 5, 6, 7), it was possible to approximate TTsh and viau with 
reasonable agreement with the in situ data. 


5. Conclusions 

We introduced the blast wave regression (BWR) technique. This technique re¬ 
constructs the evolution of fast ICMEs/shocks associated with type H radio 
bursts. The technique assumes that interplanetary shocks propagate as blast 
waves through a stationary homogeneous ambient solar wind, whose density 
decays with the square of the heliocentric distance. Using frequency data from 
the TII emission, the technique develops a dispersion relation whose regression 
solves the kinematics (speed and position) of the IP shock. We applied the 
BWR technique to analyze 8 ICME/shock events associated with TII emissions 
occurred during the decaying phase of solar cycle 23. The results suggest that 
the shock events propagated as blast waves beyond ^21 Rq. The BWR results 
also showed notable consistencies with data from coronagraph images, speeds 
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deduced from IPS observations and in situ measurements. The BWR calculated 
transit times and arrival speeds (at 1 AU) of shocks showed about 90% of 
consistency with the in situ data. The accuracy of this technique depends on 
the availability and quality of the type II spectral data associated with the 
ICME/shock. The results may be mislead by contamination of radio data and 
directivity issues. However, even for those unfavorable conditions, the technique 
was capable to estimate transit times and arrival speeds that differed around 
20% from their in situ counterparts. The agreement that our results showed 
with different data sets, preliminarily indicates that the BWR technique could 
be a useful tool for ICME/shock tracking and arrival forecasting, although the 
performance analysis of the latter is not under the scope of this paper. 


Appendix 

Uncertainties 


We can estimate the uncertainties associated with the BWR technique by using 
the differentiation method. Then, the uncertainty (Svi) of is given by: 


Svi _ SniAu ^ 
Ti niAu /■ 


( 8 ) 


where 6niAu and are the uncertainties of uiau and /^, respectively. Here 
we can estimate Sf^ ^ AfijA. Finally, the uncertainties on shock trajectory are 
expressed by: 
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and 


( 9 ) 

( 10 ) 

( 11 ) 


Where 5wsw , 5t and 5vq are the errors of solar wind speed, r and initial speed of 
shock front, respectively. Additionally, 5niAu is the uncertainty related to uiau 
and Sro the uncertainty of rg. 
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